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Abstract 

We demonstrate that the collinear expansion fails in the case of gluon emission from a 
fast quark produced in eA DIS. In this approximation the N = l rescattering contribution 
to the gluon spectrum vanishes. We show that the higher-twist approach by Guo, Wang 
and Zhang [4, 5] based on the collinear expansion is wrong. The nonzero gluon spectrum 
obtained in [4, 5] is a consequence of unjustified neglecting some important terms in the 
j> ■ collinear expansion. 

<n: 

1. Since the early nineties considerable theoretical efforts have been made to study the 
induced gluon emission from fast partons due to multiple scattering in cold nuclear matter 
and hot quark-gluon plasma (QGP). There have been developed several approaches to this 
phenomenon. The well known BDMPS [1] and GLV [2] approaches are based on the time- 
ordered diagram technique. The BDMPS formalism is valid for massless partons in the 
j> ! limit of strong Landau-Pomeranchuk-Migdal (LPM) [3] suppression when the number 

of rescatterings N ^$> 1. The GLV approach accounts for small number of rescatterings 
(N < 3) and applies only to thin plasmas. The higher-twist method due to Guo, Wang and 
Zhang (GWZ) [4, 5] is based on the Feynman diagram formalism and collinear expansion. 
It includes only N = 1 rescattering. The formalism has originally been derived for the 
gluon emission from a fast quark produced in eA DIS. The light-cone path integral (LCPI) 
formalism [6] (see also [7, 8, 9]) is free from the restrictions of the approaches [1, 2, 4, 5]. 
It accurately treats the mass effects, and applies at arbitrary LPM suppression. 

The LCPI [6] and BDMPS [1] approaches become equivalent at strong LPM sup- 
pression [8, 9]. The predictions of the GLV [2] formalism can be obtained in the LCPI 
approach [6] by expanding the gluon spectrum in the density of the medium. However, 
the status of the GWZ approach [4, 5] is not clear. The BDMPS-GLV-LCPI approaches 
in their original formulations neglect the quantum nonlocality in production of the fast 
partons. In the case of gluon emission in the QGP produced in AA-collisions this approx- 
imation is justified by the fact that the quantum nonlocality is ~ 1/E (E is the energy 
of the fast parton), and can be safely neglected at large E. In the GWZ approach to eA 
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DIS the nonlocal fast quark production and gluon emission have been treated on even 
footing. However, as we will show below in the applicability region of the GWZ formalism 
the quantum nonlocality in the quark production in eA DIS, similarly to AA-collisions, 
is not important for gluon emission. For this reason one could expect that the GWZ 
gluon spectrum should coincide with N — l contribution in the LCPI formalism. But the 
two spectra are not identical. Indeed, the GWZ gluon spectrum for the Gaussian nuclear 
density n^r) oc exp (—r 2 /2R\) at z <C 1 (hereafter z = uj/E, where uj is the gluon energy 
and E is the struck quark energy) can be written as (up to unimportant factors) 
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where Pc q = Cj?[l + (1 — z) 2 ]/z, i< 1, Gn is the nucleon gluon density. To leading order 
in a s xGn{x,Q 2 ) ~ 3oi ^f \ n [Q 2 / y 2 ) (Q 2 ~ p 2 ^ ^ [ s infrared cutoff). On the other hand 
for the N — l gluon spectrum in the LCPI approach one can obtain for massless partons 
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One can see that the integrand in (2) does not have any logarithmic factor which could 
be interpreted as the nucleon gluon density entering (1). 

In the present paper we clarify the situation with the discrepancy between the predic- 
tions of the GWZ and LCPI approaches. We demonstrate that the approximations used 
in [4, 5] really lead to a disagreement with the LCPI approach. However, contrary to the 
results of [4, 5] the consistent use of the method of [4, 5] gives a vanishing gluon spectrum. 
This fact is a consequence of failure of the collinear expansion. We show that the authors 
of Refs. [4, 5] obtained the nonzero spectrum just due to unjustified neglecting some 
important terms in the collinear expansion. 

2. As in [4, 5], we consider eA DIS (we will discuss the case when the virtual photon 
strikes out a quark with charge e q ). As usual we choose the virtual photon momentum in 
the negative z direction, and describe the and 3 components of the four- vectors in terms 
of the light-cone variables y ± = (y° ±y 3 )/v2. In the laboratory frame the photon energy 
reads v = Q 2 j2m^XB {%b is the Bjorken variable). The p~ momentum of the struck quark 
equals y / 2z/(l + Xs/2mArz/). The transverse momentum integrated distribution for the gq 
final state in eA DIS can be described in terms of the semi-inclusive nuclear hadronic 
tensor dW^ /dz. To leading order in v the spin effects in the final-state rescatterings of 
fast partons can be neglected. It ensures that the spin structure of dW^ /dz is the same 
as for the usual hadronic tensor W^ in eN DIS. It allows one to write the semi-inclusive 
nuclear hadronic tensor as dW^ 1 ' /dz = —e^g^dfA/dz, where dJA/dz is the semi-inclusive 
quark distribution of the target nucleus. Neglecting the EMC and shadowing effects it 
can be written as (we suppress the arguments xg, Q and z) 

where df^{v)/dz is the in- medium semi-inclusive quark distribution for a nucleon located 
at r, and n A (r) is the nucleus number density. 



Let us first discuss the evaluation of the in-medium semi-inclusive quark distribution 
in the LCPI formalism. We will only concentrate on the aspects which are important 
for comparison with the higher-twist method. Particularly interesting is the question on 
the quantum nonlocality of the fast quark production in eA DIS. Our consideration will 
be physical and diagrammatic. A more formal analysis will be given in further detailed 
publication. In the LCPI approach [6] the matrix element of the q — > gq' in-medium 
transition is written in terms of the wave functions of the initial quark and final quark 
and gluon in the nucleus color field (we omit the color factors and indices) 

(gq'\S\q) = tg J dy^y^A^^y) . (4) 

Each quark wave function in (4) is written as ip(y) = exp(—ip~y + )u\(fi(y~, y T )/^2p~, 
where A is quark helicity, u\ is the Dirac spinor operator. The y~ dependence of the 
transverse wave functions is governed by the two-dimensional Schrodinger equation 



dy~ l 2/2 

with the Schrodinger "mass" /j, — p~. The wave function of the emitted gluon can be 
represented in a similar way. 

The y~ evolution of the transverse wave functions can be written in terms of the 
Green's function for the Schrodinger equation (5) 

0(2/2", j/r,2) = djk,ifc(ifr2,y2\ikA,yT)<i>(yT,yT,i)- ( 6 ) 



It allows one to represent the matrix element in terms of the transverse Green's functions 
/C. In the LCPI approach [6] the path integral representation for the Green's functions is 
a crucial step in the final stage of calculations. It allows to perform averaging over the 
medium potential at the level of the integrand in the double path integral representation 
of the gluon spectrum. This method turns out to be very convenient in accounting for 
arbitrary number of rescatterings. But for comparison with the higher-twist method we 
need only the N — 1 rescattering contribution. It can be analyzed without using the 
path integral representation of the Green's functions. For calculation of the N = 1 term 
it is enough to have the perturbative expansion of the Green's functions to the second 
order in the external potential. One can show that for gauges with potential vanishing at 
large distances (say, covariant gauges, or Coulomb gauge) one can ignore the transverse 
component At- Then, the second order expansion of K, can be written as 

fc(yT,2,y2\yT,i,Vi) = K (yT,2,y2~\yT,i,yT) -ig J dy~dy T K(y Tt2 ,y2\yT,y~)A + (y T ,y~) 

xK(y T ,y~\yT,i,yT) ~ 9 2 J dy~dy T J dz~dz T K(y Tt2 ,y2\yT,y~)A + (y T ,y~) 

y->z- 

xK(y T , y"\z T , z~)A + (z T , z~)K(z T , z~\y T ,u Vi) , (7) 

where K is the free Green's function when the external field is absent. The expansion (7) 
allows one to represent diagrammatically the N = l induced gluon emission in eA DIS by 



the set of diagrams like shown in Fig. 1 in which the horizontal solid line corresponds to 
K (— >) and K* (<—), the gluon line shows the correlator (A + (yi)A + (y 2 )) in the nucleus, 
and the vertical dashed line shows the transverse density matrices of the final partons at 
very large y~. The graphs like Fig. la describe the real induced q — > gq transition, and 
the virtual graphs like Fig. lb, coming from the interference of the second and zeroth 
order terms on the right-hand side of (7), correspond to the unitarity correction due to 
interference of of the induced and vacuum gluon emission. 

The typical difference in the coordinate y~ for the upper and lower 7*gg vertices 
in Fig. 1 (which gives the scale of the quantum nonlocality of the fast quark pro- 
duction) is given by the well known Ioffe length Lj = 1/ijinXb- For the usual nu- 
cleon quark distribution Lj is the dominating scale in the Collins-Soper formula [10] 
f N = _L / dy~e iXBP+v ~ {N\ip(-y- /2)7+-0(y~ /2)\N). Contrary to the final state with sin- 
gle quark for the gq final state the integration over the y~ coordinate of the ^*qq vertex 
is now affected by the integration over the positions of rescatterings and the q — > gq 
splitting. But for moderate xb when Lj <C Ra one can neglect the effect of rescatterings 
on the integration over y~ . One can easily show that for production of the final gq states 
with M 2 <C Q 2 the restriction on y~ from the splitting point can also be ignored. Indeed, 
the typical scale in integrating over the splitting points is given by the gluon formation 
length Lf ~ vjM 2 which is much bigger than Li at M 2 <C Q 2 (this is true for both the 
vacuum DGLAP and the induced gluon emission). Also, at Lj <C Lf one can take for 
the lower limit of the integration over the splitting points for the upper and lower parts 
of the diagrams in Fig. 1 the position of the struck nucleon. After this simplifications the 
quark production and gluon emission become mutually independent, and the df^i^/dz 
can be approximated by the factorized form 

df N (v) _ dP(r) 

JN- 



dz dz 

where the quark distribution /jv stems from the left parts of the diagrams, and the gluon 
spectrum dP/dz is described by the right parts of the diagrams evaluated neglecting the 
quantum nonlocality of the fast quark production. 

One can make some more simplifications. If one neglects the multiquark configurations 
in the nucleus both the vector potentials in the gluon correlators should belong to the 
same nucleon. For this reason the typical separation of the arguments in the correlators 
is of the order of the nucleon radius, Rn- The restriction |yf — y 2 \ < Rn allows one to 
replace the fast parton propagators between the gluon fields in the graphs like Fig. lb by 
5 functions in impact parameter space. This approximation is valid for parton energy ^> 
1/ Rn- It follows from the Schrodinger diffusion relation for the parton transverse motion 
p 2 ~ L/E. Also, the smallness of the fast parton diffusion radius at the longitudinal 
scale ~ Rn allows one to replace in other transverse Green's functions the y~ coordinates 
by the mean values of the arguments of the vector potentials in the gluon correlators. 
This approximation corresponds to a picture with rescatterings of fast partons on zero 
thickness scattering centers (nucleons). The inequality uo ^> 1/ Rn for the emitted gluon 
is equivalent to Rn <C Lf. For this reason, in the picture of thin nucleons the contribution 
of the graphs like Fig. lc,d with gluon correlators connecting the initial quark and final 
quark or gluon can be neglected since they are suppressed by the small factor Rn/ Lf. 



These approximations have been used in the original formulation of the LCPI approach 
[6] (the BDMPS [1] and GLV [2] approaches use them as well). Note that (similarly to 
the case of QGP [11]) each gluon correlator appears only in the form of an integral over 
Ay~ = y 2 — y± and at y 2 — yt- One can easily show that this ensures gauge invariance 
of the result (to leading order in a s ). 

The starting point of the higher-twist approach [4, 5] to the gluon emission in eA DIS 
is the representation of the nuclear hadronic tensor for the 7*A — ► gqX transition in terms 
of the diagrams like shown in Fig. 2. The lower soft part is expressed in terms of the 
matrix element (A\ip(0)A + (yi)A + (y 2 )i/j(y)\A), and the upper hard parts are calculated 
perturbatively. In the calculation of the hard parts to leading order in the struck quark 
energy the integrations over the y + coordinates of the vertices give conservations of the 
large p~ momentum components of fast partons. At the same time the integration over 
the p~ momentum components of the final partons ensures that to leading order in the 
struck quark energy all the y + coordinates in the soft part can be set to zero. Due to 
conservation of the large p~~ momenta of fast partons in the Feynman propagators only the 
Fourier components with p~ > are important. It means that the Feynman propagators 
are effectively reduced to the retarded (in y~ coordinate) ones. 

One can easily demonstrate that the Feynman diagram treatment of [4, 5] is equivalent 
to that in terms of the transverse Green's functions. Indeed, using the representation 

v( ~ -, - -n • f dp + dp T exp [-ip + (y 2 - y-i) + ipr(yT,2 - #r,i)] {n , 

K(vr»v 2 \vr, 1 ,y 1 )=iJ -^- ^ y | (9) 

one can write the retarded quark propagator as 
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Here u\ and u\ act on the variables with indices 2 and 1, respectively. The last term in 
(10) is the so-called contact term. It does not propagate in y~ and can be omitted in 
calculating the nuclear final-state interaction effects for fast partons. One can easily verify 
that using (10) and a similar representation for the gluon propagator the hard parts in 
the higher-twist method can be represented in terms of the transverse Green's functions 
as it is done in the LCPI treatment. Thus, the hard parts in the approach [4, 5] should 
agree with that of the LCPI formalism. 

3. The calculation of the diagrams like shown in Fig. 1 is simplified by noting that the 
free transverse Green's function for a parton with mass m can be written as 

K(y T ,2, y 2 \yT,i, i/f) = %2~ - Vi) J2 $vt{vt,i, y 2 )<t>* PT (yT,u yT) , ( n ) 

PT 

where 4> PT (yT,y~) = exp [ipryr — iy~{Pr + m 2 )/2p~} is the plane wave solution to the 
Schrodinger equation for A^ = with the transverse momentum pr- It allows one 
to represent the upper and lower parts of the diagrams shown in Fig. 1 in the form 



/ dy~dyT(fiqi(yT, y~)<f>l(yT, y~)(pq(yT, y~) where the outgoing and incoming wave functions 
have the form of the plane waves with sharp change of the transverse momentum at the 
points of interactions with the external gluon fields. This method has previously been 
used in [12] for investigation of the kinematical effects. If one ignores the finite kinemati- 
cal limits it reproduces the first order in the density term of the full LCPI expression for 
the induced gluon spectrum obtained in [13]. We have checked that all the hard parts 
evaluated with the help of the plane waves agree with that obtained in Refs. [4, 5] . 

The sum of the complete set of the diagrams contributing to the N — l spectrum can 
be written as [13, 12] 
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Here da(z,^)/dz is the cross section of gluon emission from a fast quark produced at 
distance £ from the scattering nucleon. At z <C 1 (we consider the soft gluon emission 
just to simplify the formulas) for massless partons it reads [13, 12] 
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(14) 
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Here the limit x — > is implicit, dG(k^, x)/dkr is the unintegrated nucleon gluon density 
which in leading order in a s at x <^ 1 reads 
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dpdy T exp {-ik T p){m N \V yT W a {y T + p)V yT W a (y T )\^ N ) . (15) 



Here W a (yr) = J dy~ A +a (y~ , y T ) (the color index, a, is shown explicitly), and ^n is 
the internal nucleon wave function normalized to unity. One can show that the formula 
(15) being rewritten through the nucleon wave functions with relativistic normalization 
(N'\N) = 2p + (2n) 3 5(p + —p + )5{p T —pr) is reduced to the Collins-Soper definition [10] of 
the gluon density. Eq. (15) can also be written as 

dG(k 2 T ,x) N 2 C 
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where a(p) is the well known dipole cross section given by 
Sna 2 s CF 
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Using (12)-(15) with the phase W a calculated in the Born approximation one can obtain 

(2). 

The collinear expansion corresponds to replacement of the hard part by its second 
order expansion in kx (we suppress all the arguments except for kr) 
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Then, to logarithmic accuracy da(z,^)/dz oc J dprxG(pj',x)Vl T H(pT, kr, z, Olj£ =tr But 
from (14) one can easily obtain V^ H\% =0 = Q It is also seen from averaging of the hard 
part over the azimuthal angle of kr which gives explicitly 

1 r -> 6{hr — Pt) 
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Thus, contrary to the expected in the collinear approximation dominance of the region 
kr 1z> Pt only the region kr > Pt contributes to the gluon emission, and formal use of the 
collinear expansion gives completely wrong result with zero gluon spectrum. 

The fact that the N — \ gluon spectrum vanishes in the collinear approximation agrees 
with prediction of the harmonic oscillator approximation in the LCPI [6] and BDMPS 
[1] approaches. The full gluon spectrum in these approaches is expressed in terms of the 
Green's function of the Schrodinger equation with an imaginary potential which is pro- 
portional to the dipole cross section. The corresponding Hamiltonian takes the harmonic 
oscillator form for a quadratic parametrization a(p) = C p 2 . The N — 1 contribution 
to the gluon spectrum in the oscillator approximation should coincide with prediction of 
the collinear expansion in the GWZ treatment [4, 5] since the oscillator approximation in 
the LCPI and BDMPS approaches is equivalent to the collinear approximation of [4, 5]. 
Indeed, the quadratic form of the dipole cross section corresponds to the vector potential 
approximated by the linear expansion A + (y~,y T + p) m A + (y~,y T ) + pV yT A + (y~,y T ) 
which can be traced back to the collinear expansion in momentum space. For a target of 
thickness L in the oscillator approximation the BDMPS and LCPI approaches for massless 
partons give the spectrum [1, 8, 13] 

dP OA _ a s P Gq {-" 



dz n 



ln|cos£lL|, (20) 

where f2 = J —iC^n/ z{l — z)E, C3 = CCa/Cf- Keeping only the first term in the 
expansion of the spectrum (20) in the density one gets [13] 
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Since the right-hand side of (21) oc n 2 it corresponds to the N = 2 rescatterings [13], and 
the contribution of N — 1 rescattering is absent. Note that from the point of view of the 
representation (13) the zero N — 1 gluon spectrum in the oscillator approximation is a 
consequence of the fact that in this case dG/dkx oc <5(A>r) (as one sees from (16)). 

4. Let us now discuss why the calculations of Refs. [4, 5] give nonzero gluon spectrum. 
In [4, 5] the nonvanishing second derivative of the hard part comes from the graph shown 
in Fig. 2b (at z <^i 1). The authors use for the integration variable in the hard part of 
this graph the transverse momentum of the final gluon, It- The /r-integrated hard part 
obtained in [5] (Eq. 15 of [5]) reads (up to an unimportant factor) 

H(k T ) oc / T R(y~, </r, V2, h, k T ) , (22) 
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is an analog of the last factor in the square brackets in (14) (our z corresponds to 1 — z 
in [4, 5]), the coordinates y~, y± 2 correspond to the quark interactions with the virtual 
photon and t-channel gluons. In calculating V\ T H the authors differentiate only the 
singular factor 1/(It ~ kr) 2 - But the omitted terms from differentiating the factor R are 
important. After the It integration they almost completely cancel the contribution from 
the l/ilr — kT) 2 term. Indeed, after putting yf = y^ and changing the integration variable 
It — » (It + kr) the right-hand part of (22) does not depend on kr at all. If one does not 
put yf = y 2 i there will be some nonzero contribution to Wl T H\^ =Q which, however, is 
suppressed by the small factor (Rn/ ' Lf) 2 . Keeping such contributions does not make any 
sense since they are clearly beyond predictive accuracy of the approximations used in 
[4,5p. 

5. In summary, we have demonstrated that the collinear expansion fails in the case of 
gluon emission from a fast massless quark produced in eA DIS. It gives a zero N — l rescat- 
tering contribution to the gluon spectrum. This agrees with vanishing N — l contribution 
to the gluon spectrum in the BDMPS [1] and LCPI [6] approaches in the harmonic oscilla- 
tor approximation which corresponds to the collinear expansion in momentum space. The 
nonzero gluon spectrum obtained in [4, 5] is a consequence of unjustified neglecting some 
important terms in the collinear expansion. The established facts demonstrate that the 
GWZ approach [4, 5] is completely wrong. Its predictions for eA DIS and jet quenching 
in AA collisions do not make sense. 
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